Study of anomalous mobility of polar molecular solutions 
by means of the site-site memory equation formalism 
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In this work, the memory equation approach is applied for theoretical study of dynamics of 
polar molecular liquids described by the interaction site model. The study includes the tempe- 
rature-density(pressure) dependence of the translational diffusion coefficients D and orientational 
relaxation times r for infinitely dilute solutions of acetonitrile and methanol in water, and methanol 
in acetonitrile. Calculations are performed over the range of temperatures and densities employing 
the SPC/E model for water and optimized site-site potentials for acetonitrile and methanol. Despite 
an approximate character of the model potentials and closure relation used, the theory is able 
to reproduce qualitatively all main features of temperature and density dependences of D and r 
observed in computer and real experiments. In particular, anomalous behavior, i.e. the increase in 
mobility with density (pressure), is obtained for D and r of methanol in water, while acetonitrile in 
water or methanol in acetonitrile do not show deviations from the usual. The observed enhancement 
in the molecular mobility is interpreted in accordance with the concept by Yamaguchi et al. [J. Chem. 
Phys. 119 (2003) 1021], i.e. in terms of two competing origins of friction, which interplay with each 
other as density increases: the collisional and dielectric frictions that have tendency, respectively, 
to strengthen and weaken with increasing density. 
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I. INTRODUCTION 

The enhancement of the molecular mobility of polar 
molecules in dense dipolar liquids under pressure, called 
anomalous molecular mobility, has recently attracted 
considerable attention in solution chemistry experiments 

and simulations HEHSHH0I1I!. Conventional 
continuum-based theories of liquid matter collapsed when 
applied for the description of anomalous molecular be- 
havior, so that to make theoretical justification of the 
observed processes is topical from the point of view sta- 
tioned on the statistical mechanical approach |lC|. In 
order to clarify the microscopic structure of water, which 
is a typical representative of a liquid with anomalous be- 
havior, and its mobility, it has been studied recently 
by Yamaguchi et al. using the memory equation ap- 
proach for molecular liquids described by the interaction 
site model The authors concluded, that anomalous 
pressure dependence of dynamic properties of water can 
be explained by the strong interatomic Coulomb interac- 
tion correlated with density fluctuations. In particular, 
when the number-density fluctuation in the low wave- 
number region decreases with increasing pressure, the 
Coulomb coupling between the low wave-number charge- 
density fluctuation and hydrodynamic dielectric mode is 
reduced. In the case of another polar liquids, like acetoni- 
trile, the effect of Coulomb interaction can be blocked by 
even more substantial repulsive interaction, which makes 
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it dissimilar with water. Phrasing it differently - anoma- 
lous molecular mobility strongly depends on the molec- 
ular geometry and dominant type of interatomic inter- 
action, and more studies is required for better under- 
standing of its nature. In this paper we pay attention to 
the study of anomalous pressure behavior of polar binary 
molecular solutions by considering such popular models 
as methanol in water, acetonitrile in water and methanol 
in acetonitrile. In particular, we are interested in the den- 
sity (pressure) dependence of translational diffusion coef- 
ficients and reorientation relaxation times of solutes. The 
rest of the paper is organized as follows. In section ITT1 we 
write down equations of motion for site-site intermediate 
scattering functions, and give relations for translational 
diffusion coefficient and reorientation correlation time in 
terms of these functions. Equilibrium properties of the 
system are obtain using DRISM, and memory kernels are 
estimated by the mode-coupling approximation. Set up 
for systems under consideration and details of numerical 
procedures are described in section ITTll Obtained results 
are discussed in section llVl and conclusions are given in 
section Ivl 



II. THEORY 

For the description of time-dependent properties of the 
system, we will follow the formalism by Hirata and Chong 
12], which is the unification of the theory of dynamical 
processes in simple liquids based on the Mori approach 
|l3l fl4| and reference interaction site model for molec- 
ular liquids [TEl ITtl lri| . In their formalism, parameters 
of shortened description of the system are partial num- 
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ber densities and longitudinal current densities. Then, 
of practical interest are elements of the matrix of the 
site-site intermediate scattering functions F(k;t) and the 
matrix of their self-parts F s (k;t) defined by 

F^{k-t) = i<p».*(k;(V(k;i)>, (la) 
*T(M) = l<p a >*(k;0Mk;i)) S) (lb) 

where N is the total number of particles, p a (k;t) is the 
site a number density in reciprocal space, k is the wave- 
vector, k = |k|, t is time, * means complex conjuga- 
tion, angular brackets denote appropriate statistical ave- 
rage (e.g., canonical Gibbs ensemble average), and suffix 
"s" stands for "self" and means correlations between two 
sites of a same molecule. In the limit of infinite dilution, 
which is considered here, the solution density is essen- 
tially determined by the solvent. Consecutively, one has 
two types of memory equations for the solvent subsys- 
tem and one for the solute subsystem (indicated by the 
superscript "u") which are: 



F(fc;t) 



{u>t)?{k-t)- / drK(fc;r)F(fc;t-r), 
Jo 



(2a) 



F B (fc;t) = -<«k>.F.(M)- / dr K s (fc;r)F s (M-r), 

Jo 

(2b) 

(u; 2 ) b u F s u (M) ~ / drK s u (fc;r)F s U (fc;t-r). 
Jo 

(2c) 



F,(M) 



In these equations, dot over the quantity means its time 
derivative; the memory function matrices, denoted as 
K(fc;i), K B (k;t) and K" (fc;t), describe the friction on the 
motion of interaction sites; quantities (u> 2 ), (w|) s and 
(w?)" are normalized second order frequency matrices 
given by 



,2\u 



fc 2 J(fc)S _1 (fc), 
fc 2 J(fc)S s - x (fc), 
k 2 r(k)S«>-\k), 



(3a) 
(3b) 
(3c) 



where S{k) = F(k;t = 0), S s (fc) = F s (k;t = 0) and 
S" (k) = Fg (k; t = 0) are matrices of static site-site struc- 
ture factors and their self parts, respectively, while J (A;) is 
the matrix of static longitudinal site-current correlation 
functions 



i,3 



— 2k-(r" — r7 ) 



(4) 



where subscripts i, j refer to molecules, rf = vf (0) is the 
initial position of site a, and vf z = vf z (0) is longitudinal 
component of the initial velocity of site a. The analyti- 
cal expression of J (k) for arbitrary shape of the molecule 



has been given recently by Yamaguchi et al. [lRj . Defini- 
tion for J u (fc) is similar to the one given by equation 
with the difference that summation runs over the solute 
molecule only. Initial values of intermediate scattering 
functions, that we need to solve memory equations (0, 
can be obtained using RISM theory [H El E3- That 
theory predicts static structure of molecular fluids via the 
calculation of their site-site pair correlation functions. In 
this work, however, in order to have a better descrip- 
tion for the dielectric susceptibility of solvent we employ 
the formalism of DRISM that uses experimental value of 
dielectric constant 0, ■ 

Quantities of our interest in this paper are solute's 
translational diffusion coefficient D and rank-1 reorien- 
tation correlation time r. Following the derivation pro- 
cedure presented in [2lJ , the translational diffusion coef- 
ficient D is obtained as 

1 



D = - 



dt Z a '<(t) 



3./o 
— lim 



(5) 



where Z a ~<(t) is the site-site velocity autocorrelation 
function with sites a and 7 belonging to the same 
molecule. The rank-1 reorientation correlation function 
C M (£) is defined by 



C M (t) = 



(/*«(Q)/ii(*)) 



(6) 



where Hi(t) is a vector fixed on the molecule i. In our 
case it is the dipole moment and therefore can be de- 
scribed by the linear combination of site coordinates as 
MiW = Ea z a ra W wrtn z a being site partial charges. 
Putting that into equation JSJ and using properties of 
time-correlation functions 22] one arrives at 



£i£q 7 *aMr, ? (0)r7(*)) 



(7a) 
(7b) 



Hence, the time development of both Z a ^(t) and C M (i) is 
governed by the memory equation for the self-part of the 
site-site intermediate scattering function F"(fc;i). The 
rank-1 reorientation relaxation time is defined as 231 



dtC„(t). 



(8) 



Since the dynamics we are interested in corresponds to 
the long-time limit, memory kernels for memory equa- 
tions can be constructed using the mode-coupling ap- 
proximation [23j . In works by Chong et al. |2J, |25j the 
conventional mode-coupling theory has been extended to 
the case of molecular liquids based on the interaction-site 
model. It has been shown, however, that the proposed 
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expressions for memory functions underestimate friction 
in orientational motions |26(. In our study we use the 
recipe by Yamaguchi and Hirata |2g, who suggested a 
modified expression that includes the interaxial coupling. 



III. DETAILS OF THE MODEL SETUP AND 
NUMERICAL PROCEDURES 

We performed calculations for acetonitrile in water, 
methanol in water and methanol in acetonitrile, all in 
the case of infinite dilution. As for the structure and the 
intermolecular potential of water we employed a model 
of the extended simple point charge (SPC/E) [2?|. We 
also put the Lennard-Jones (LJ) core on the hydrogen 
atoms in order to avoid the undesired divergence of the 
solution of the DRISM integral equation. The L J param- 
eters of the hydrogen atom, the depth of the well and the 
diameter, were chosen to be 0.046 kcal/mol and 0.7 A, 
respectively. 

In acetonitrile and methanol the methyl group was 
considered to be a single interaction site located on the 
methyl carbon atom. So that all chemical compounds 
consisted of three sites which interact through the pair 
potential [HHU 



( r ij ) — 4e Q - 
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u a7 



(9) 



i.e., LJ plus Coulomb. Here r^ 7 = \rf — rj|; parame- 
ters e Q7 and a a j are LJ well-depths and LJ diameters 
defined as e Q7 = ^/e Q e 7 and cr Q7 = (a a + c 7 )/2, respec- 
tively, with <j a being the LJ diameter of a single site. 
Point charges for acetonitrile were chosen to reproduce 
electrostatic potential obtained in ab initio calculations 

In calculations for acetonitrile or methanol in water 
temperature of the system was varied from 258.15 to 
373.15 K, and density of water from 0.9 to 1.2 g/cm 3 ; for 
the case of methanol in acetonitrile temperature of the 
system was varied from 293.15 to 323.15 K, and density 
of acetonitrile from 0.6726 to 0.815 g/cm 3 . Connection of 
the water parameters with thermodynamic pressure can 
be established, e.g., using the multi-parametric equation 
of state [30J (except for the metastable regions where re- 
liable data is lacking). 

Temperature/density dependent dielectric constant e 
for water used in numerical calculations has been eval- 
uated as a physical solution of an empirical nonlinear 
equation [31j: 



1 
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9.32 • 10 4 (1 + 



153 ' 



(1 - 3/ufT 



(10) 

where v is a molar volume in units of cm 3 /mol, and T 
is thermodynamic temperature in K. This equation has 
been used also in such temperature/density points where 



no experimental values exist. For density and dielectric 
constant of acetonitrile at different temperatures we used 
experimental values indicated in Table [I] 



TABLE I: Density and dielectric constant for acetonitrile as 
functions of temperature: experimental data used in our com- 
putation. Temperature is in K, and density is in g/cm 3 . 



293.13 295.05 



.782 a 
38. 8 C 



37. 5 d 



298.15 
.7762 b 
36.69 c 



303.15 308.15 313.15 323.15 



.7712 a 
35.93 c 



.7652° 



.7603 a 



.7492 a 
33.5 f 



a l3l. b l3l. C [3H. d t35l. e l3l. f Extrapolated 

In the numerical part of this study we calculated first 
the site-site static structure factor by solving the DRISM 
equation using the intermolecular interaction, molecular 
shape, temperature and density. In order to improve 
the convergence of the DRISM calculation, we used the 
method of the modified direct inversion in an iterative 
space (MDIIS) proposed by Kovalenko et al. .37}. The 
value of the parameter of DRISM theory a has been set to 
0.1 A 2 . From the static site-site structure factor, we cal- 
culated the site-site intermediate scattering function us- 
ing the site-site memory equation with the mode-coupling 
approximation for the memory kernels. The memory 
equation has been time-integrated numerically. Time- 
development of correlation functions in the k —* limit 
was treated separately by the analytical limiting proce- 
dure of theoretical expressions. In the numerical pro- 
cedure, the reciprocal space was linearly discretized as 
k = (n + ^)Afc, where n is an integer from to — 1. 
Values of Ak and Nk are 0.061 A -1 and 2 9 = 512, respec- 
tively. The choice for N). as the power of two came as the 
requirement of the subroutine for the fast Fourier trans- 
form, which has been used in DRISM / MDIIS. The dif- 
fusion coefficient D was calculated from the asymptotic 
slope of the time dependence of the mean square dis- 
placement, and the orientational relaxation time r was 
obtained by time-integration of the rotational autocorre- 
lation function. 



IV. RESULTS AND DISCUSSIONS 
A. Structural properties 

When applied to the liquid state, external pressure 
induces changes both in its structure and dynamics of 
molecules. The anomalous behavior of the molecular mo- 
bility with pressure, i.e. dynamics, has been assigned to 
the hydrogen bonding properties of the system, which in 
its turn can be related to structural properties. The idea 
is that the hydrogen bond or hydrogen bonding network 
are distorted upon compression, so that the molecule has 
less hydrogen bonds at higher pressure, which makes its 
translational and rotational motion looser. On the other 
hand, extensive compression results in significant grow 
of repulsive forces between molecules. In fact, the liq- 



4 





FIG. 2: Water-methanol site-site radial distribution functions at T = 273.15 K and set of densities from 0.9 to 1.125 g/cm 3 for 
water, obtained by the DRISM/HNC integral equation theory. In the notations used, first site always belongs to water, and 
second site always belongs to methanol. Arrows show directions of alternations due to an increase in pressure. 



uid structure is determined essentially by these two com- ular constraints or molecular geometry. There was much 
peting factors: hydrogen bonding and repulsive interac- attention to this topic in scientific literature and we shall 
tion, which interfere each other through the intramolec- not repeat results of these studies here. Instead of that 
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FIG. 3: Water-acetonitrile site-site radial distribution functions at T = 273.15 K and set of densities from 0.9 to 1.125 g/cm 3 
for water, obtained by the DRISM/HNC integral equation theory. In the notations used, first site always belongs to water, and 
second site always belongs to acetonitrile. Arrows show directions of alternations due to an increase in pressure. 



we can focus our analysis on some properties related to 
the subject. Figs. ^"Elshow site-site RDF's for neat wa- 
ter, water-methanol, water-acetonitrile and acetonitrile- 
methanol infinite dilution solutions, respectively, calcu- 
lated by the DRISM/NHC theory for set of densities at 
constant temperature. Despite the comparatively big 
number of site-site combinations, there are some com- 
mon features in the behavior of each individual site-site 
RDF, that can be categorized as follows. First, one sees 
that there are RDF's with the first peak either remained 
almost unchanged, or enhanced with the increase of den- 
sity (pressure). Second, there are first peaks that are lo- 
cated at positions much shorter, than other peaks, and 
that all of them subside with density(pressure) a lot. The 
former case is a reflection of packing effect. It is because 
of the result of pressure pushes the next nearest neighbor 
toward the central molecule and makes the first coordina- 
tion shell to be outlined more clearly, which is also con- 
firmed by the drift of majority of peak locations into di- 
rection of shorter separations. The effect of pressure also 
distorts the preferable mutual orientation of molecules, 
which docs not favor hydrogen bonds. The latter case 



in given previously peak categorization is an example of 
such state of things. To be more specific, first peaks in 
Figs. Htb),Etb) andHte),Etd), andUb) all are evidences 
of hydrogen bonds in the system, and all of them subside 
with the increase of density. There is, however, a differ- 
ence in their nature. For methanol in water the hydrogen 
bonding is realized in two ways: once when the proton 
acceptor is solvent, i.e. water, Fig. |2Jb); and once when 
the proton acceptor is the solute, i.e. methanol, Fig. Ote). 
In the first case the peak amplitude is much lower, than 
in the second case, indicating that the first type of hy- 
drogen bonding may be somewhat weaker. Nevertheless, 
due to the fact that both solvent and solute can evident 
as both proton donor and acceptor, they are able to make 
either hydrogen bond network, as in the case of neat wa- 
ter, or hydrogen bond zigzag chains, as in the case of 
pure methanol or water-methanol mixture. For acetoni- 
trile in water the hydrogen bonding is realized only in one 
way: when the nitrogen site of acetonitrile serves as the 
proton acceptor, Fig. O^d). The peak amplitude of that 
hydrogen bond is also lower compared to neat water, Fig. 
H^b), or water-methanol, Fig. |2Je), but general behavior 
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FIG. 4: Acetonitrile-methanol site-site radial distribution functions at T — 293.15 K and set of densities from 0.6726 to 0.815 
g/cm 3 for acetonitrile, obtained by the DRISM/HNC theory. In the notations used, first site always belongs to acetonitrile, 
and second site always belongs to methanol. Arrows show directions of alternations due to an increase in pressure. 



of RDF with density (pressure) is essentially the same. 
The lack of ability for acetonitrile to serve as the proton 
acceptor and donor has far reaching consequences and re- 
sults in impossibility to create a hydrogen bond network 
neither in pure substance or between water an acetoni- 
trile molecules in their mixture. Thus, the solution of 
acetonitrile in water represents situation when hydrogen 



bonds are realized between both types of molecules, but 
hydrogen bond network is possible to establish between 
water molecules exclusively. Finally, for methanol in ace- 
tonitrile the only hydrogen bond available in the system 
is between the nitrogen site of the solvent and hydrogen 
site of the solute, Fig. Efb). Any type of the hydrogen 
bond network does not appear in such solution at all. 
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B. Dynamical properties 

Fig. [5] shows the density(pressure) dependence of the 
normalized translational diffusion coefficient D/Dq and 
normalized reorientation relaxation time t/tq for inves- 
tigated solutes at various temperatures. Here D and 
To are translational diffusion coefficient and reorienta- 
tion relaxation time, respectively, of the solute at the 
ambient density for the solvent, which is p — 0.997047 
g/cm 3 for water and p = 0.782 g/cm 3 for acetonitrilc. In 
particular, Fig. 0Ja) shows that for sufficiently low tem- 
perature the diffusion coefficient of methanol in water 
first increases with density (pressure), and then smoothly 
changes its behavior to the normal one, i.e. decreases 
with density (pressure). Relatively flat maximum is ob- 
served for the density slightly smaller than that at the 
ambient condition. At higher temperatures maximum is 
not observed. Similar features can be mentioned about 
behavior of the diffusion coefficient of acetonitrile in wa- 
ter, Fig. EJb). But in this case the increase of the 
diffusion coefficient at the lowest temperature is very 
small and can be monitored only by comparing num- 
bers for neighboring density points. The cubic spline 
curve, which connects these points, makes this very flat 
maximum observed visually. Behavior of the diffusion 
coefficient at higher temperatures does not deviate from 
the normal one, as in the case of methanol in water. Ex- 
perimental measurements of diffusion coefficients for ace- 
tonitrile in water and acetonitrile-D3 (CD3CN) in water, 
both at T — 303 K, reported by Nakahara et al. Q 
testify the same tendency in their density(pressure) be- 
havior as obtained from the theory. Figure |S{c) demon- 
strates monotonous decrease of the diffusion coefficient of 
methanol in acetonitrile in the entire range of densities 
and at all investigated temperatures. Fig. |£fd) is a typ- 
ical example of anomalous density (pressure) dependence 
of the reorientation relaxation time. It is clear to see that 
it first decreases with density(pressure) and then starts 
to behave normally, i.e. increases with density (pressure). 
The order of decrease smoothly varies from about 40% 
at T = 258.15 K to about 3% at T = 298.15 K. At 
T = 373.15 K anomaly is not observed, but the increase is 
very slow for the quite wide range of densities (pressures). 
In the case of acetonitrile in water, Fig. [5fe), normal be- 
havior of the reorientation relaxation time is observed 
almost entirely except for the lowest temperature, where 
t has very shallow minimum. Finally, Fig. EJf) demon- 
strates monotonous increase of the reorientation relax- 
ation time of methanol in acetonitrile in the entire range 
of densities and at all investigated temperatures. In such 
a way, it can be regarded as a typical example of nor- 
mal density (pressure) behavior of r. Filled triangles and 
filled squares in Fig. 03 are results of MD simulation fjj 
at T = 298 K for methanol in water and acetonitrile in 
water, respectively, both at the infinite dilution. And 
asterisk symbols are experimental measurements Q at 
T = 298.2 K. One can see rather satisfactory correlation 
between the results of our computation, MD simulations 



and experiment, both qualitatively and quantitatively. 
It should be noted, however, that in terms of absolute 
values of D and r agreement with simulation and exper- 
imental data is poor. 

Theoretical aspects of anomalous molecular mobility 
have been considered previously by the example of neat 
water Within the model employed, the authors pro- 
posed that the enhancement of the molecular mobility 
by compression can be related to the suppression of the 
number-density fluctuations in the low-fc region (fc is the 
wave- number), rather than only the breakdown of the 
tetrahedral hydrogen-bonding network structure of wa- 
ter. The reasoning of this statement coherently follows 
from the examination of the behavior of memory func- 
tions and is in harmony with the heuristic explana- 
tion of the picture 13 • Let us imagine a molecule rotat- 
ing in a polar media. The rotational motion will induce 
a relaxation process of surrounding molecules in order 
to make themselves aligned to the electric field produced 
by the new orientation of the molecule in concern. The 
energy dissipation associated with the relaxation process 
is an origin of the friction on the rotational motion, or 
the dielectric friction. The larger is the charge density 
fluctuation, the greater is the dielectric friction. When 
the liquid is compressed by an external pressure, the 
mechanism of electrostatic friction on the dielectric re- 
laxation is mostly the same except that molecules in the 
solution are packed more tightly. Higher packing frac- 
tion prevents large number density fluctuations leading 
to smaller non-uniformity in the the polarization density 
and, as the result, smaller heterogeneity in the charge 
density. Then the electrostatic friction on dielectric re- 
laxation is smaller and means that the dielectric relax- 
ation time becomes shorter. Such an acceleration en- 
hances the mobility of molecules in solution through the 
dielectric friction mechanism. Spherical shape of the 
molecule, or equivalently - spherical shape for the re- 
pulsive short-range interaction, favors anomalous molec- 
ular mobility. It is because the dominant contribution to 
the rotational and therefore dielectric friction is defined 
by the type of the interaction and its range, which ex- 
plains the difference in the behavior of nearly spherical 
shape molecule of methanol and rod-like shape molecule 
of acetonitrile. For methanol the dominant contribution 
comes from the long-range Coulomb interaction, while 
for acetonitrile from the short-range repulsive one. As 
the result, methanol in water exhibits more substantial 
anomaly compared to acetonitrile in water. The anoma- 
lous behavior is attributed to the strong electrostatic 
interaction - the "hydrogen bond", among the solvent 
molecules and those between solute and solvent. The 
former causes the decrease in the dielectric relaxation 
time with pressure, while the latter induces the coupling 
between the dielectric mode of the solvent and the rota- 
tion of the solute. The anomaly is largely suppressed for 
acetonitrile in water due to enhanced significance of the 
repulsive core in the molecule. Finally, methanol in ace- 
tonitrile does not show any indication of the anomalous 
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FIG. 5: Normalized translational diffusion coefficient D/Do and reorientation relaxation time t/to for for methanol in water 
(a), (d); acetonitrile in water (b), (e); and for methaol in acetonitrile (c), (f), respectively. □, A, O - theory, ■, A - results of 
MD simulation 9] for T — 298 K, * - experimental data Q] for T = 298.2 K, lines connecting open symbols are cubic splines 
for the eye-guide. 



density dependence because of strong collisional friction 
on the collective reorientation of the solvent in this case, 
so that the dielectric relaxation becomes slower with pres- 
sure. Higher density or pressure just magnifies the effect 
of the repulsive core upon rotation and because of that - 
the collisional friction. 

V. SUMMARY 

In present paper we have calculated the density de- 
pendence of the translational diffusion coefficients and 
rank-1 reorientation relaxation times for acetonitrile and 



methanol in water, and methanol in acetonitrile at var- 
ious temperatures. Calculations have been performed 
using the site-site memory equation with the mode- 
coupling approximation for memory kernels, and the 
DRISM/HNC for static properties. For simplicity of 
theoretical and computational procedures solutions have 
been considered in the limit of infinite dilution. Calcu- 
lated quantities have been obtained to behave anoma- 
lously with density(pressure) for methanol in water. In 
particular, translational diffusion coefficient D may in- 
crease for a while with density at low temperature (super- 
cooled region), and reorientation relaxation time r may 
decrease with density to form a minimum in the vicin- 
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ity of ambient condition density for water in the quite 
wide range of temperatures. Similar computations for 
acetonitrile in water show tiny anomaly in the behavior 
of D and r only at the lowest (supercooled region) tem- 
perature, while for other regions there is no deviation 
from the usual behavior. And for methanol in acetoni- 
trile there is no deviation in the entire region of inves- 
tigated densities and temperatures. This picture is con- 
sistent with results of experimental observation and MD 
simulation, also quantitative agreement is not as good as 



qualitative. The physical origin of the anomalous density 
dependence of molecular mobility is interpreted in terms 
of two competing origins of friction, which interplay with 
each other as density increases: the collisional and dielec- 
tric frictions which, respectively, strengthen and weaken 
with increasing density. Presented results are first in sci- 
entific literature to realize anomalous molecular mobility 
of polar solute in water by means of the statistical me- 
chanical theory. 
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